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^ O^ . Abstract 

^S| . We present a quantitative confocal-microscopy study of the tran- 

^ ' sient and final microstructure of particle-stabilised emulsions formed 

OO , via demixing in a binary liquid. To this end, we have developed an 

image-analysis method that relies on structure factors obtained from 
discrete Fourier transforms of individual frames in confocal image se- 
quences. Radially averaging the squared modulus of these Fourier 
transforms before peak fitting allows extraction of dominant length 

^^ ' scales over the entire temperature range of the quench. Our proce- 

dure even yields information just after droplet nucleation, when the 
(fluorescence) contrast between the two separating phases is scarcely 
discernable in the images. We find that our emulsions are stabilised on 

^\ ' experimental time scales by interfacial particles and that they are likely 

to have bimodal droplet-size distributions. We attribute the latter to 
coalescence together with creaming being the main coarsening mecha- 
nism during the late stages of emulsification and we support this claim 
with (direct) confocal-microscopy observations. In addition, our results 
imply that the observed droplets emerge from particle-promoted nu- 
cleation, possibly followed by a free-growth regime. Finally, we argue 
that creaming strongly affects droplet growth during the early stages 
of emulsification. Future investigations could clarify the link between 
quench conditions and resulting microstructure, paving the way for 
tailor-made particle-stabilised emulsions from binary liquids. 



1 Introduction 

Particle-stabilised emulsions, in which colloidal particles rather than molec- 
ular surfactants provide stabilisation vs macroscopic phase separation, have 
received increasing attention over the past two decades [1]. The reasons for 
this growing interest have been twofold: 1) particle-stabilised emulsions are 
model arrested systems and 2) they have significant potential for applications 
in the food, personal-care, pharmaceutical, agricultural and petrochemical 
industries [H El [3] . Typically, these so-called Pickering- Ramsden emulsions 
are fabricated by direct mixing of two immiscible liquids in the presence of 
colloidal particles, e.g. via limited coalescence [4J. In this process, an excess 
of liquid-liquid interface is created, i.e. there are not enough particles to 
cover all of it. Subsequently, the emulsion coarsens, thereby reducing the 
surface-to- volume ratio of the droplets, until the area fraction of interfacial 
particles is sufficient to stabilise them. Allowing superior control over droplet 
size and variance, microfluidic techniques have also been employed to fab- 
ricate Pickering-Ramsden emulsions [5]. Recently, using partially miscible 
rather than immiscible liquids, a promising (reversible) alternative for the 
fabrication of particle-stabilised emulsions has been reported — demixing 
via nucleation and growth of droplets in a binary liquid containing colloidal 
particles [6J. 

Upon quenching a binary mixture from the single-fluid phase into the 
two-phase region of its phase diagram (figure I]) El |9] , there are several 
ways in which the liquid components can separate. If the phase diagram 
is symmetric, a deep quench through the critical point will induce spinodal 
decomposition (figure [TJa)): the liquids separate via coarsening of a bicon- 
tinuous domain pattern. A shallow, off-critical quench can take the system 
into the metastable region of the phase diagram, in between the binodal and 
spinodal lines (figured]), where demixing proceeds via nucleation and growth 
of droplets. If the dynamics of the two phases operate at very different time 
scales, demixing may proceed via viscoelastic phase separation, which can 
lead to the formation of various (transient) domain patterns |1U| . 

Domain coarsening is driven by the energy cost of the liquid-liquid in- 
terface, which is quantified by the interfacial tension 7. Interfacial particles 
provide an alternative means to reduce this (energy-expensive) interfacial 
area. This is described by 

AGd = TTr^7{l-\cose\f , (1) 

where AG^j is the free energy of detachment of a spherical particle of radius r 
and contact angle 9, as measured through the more polar phase [IJ. Even for 
the partially miscible liquids hexane and methanol, which have a relatively 
low interfacial tension 7 of order 0.3 mN-m~^ at room temperature (Tfoom)! 
equation ([1]) predicts AGd ~ 6 x 10'^ kBTroom for an 'average' colloidal 
particle, i.e. 1 /im in diameter [3 [11]. In other words, colloidal particles can 
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Figure 1: (a) Schematic phase diagram of a symmetric binary hquid a/3: c/9o"* 
is the critical volume fraction of liquid /3. Grey arrow: a deep quench from 
the single-fluid phase into the unstable region of the phase diagram, inducing 
spinodal decomposition, (b) Phase diagram for hexane-methanol [7]: 99^ is 
the single-fluid volume fraction of methanol; mole fractions were converted 
into volume fractions using linear fits to the temperature-dependent densities 
of hexane and methanol ^. Grey arrow: an off-critical quench, used here 
to form particle-stabilised emulsions via nucleation and growth of droplets. 
Panel (b) adapted by permission of The Royal Society of Chemistry [9] . 



become irreversibly attached to liquid-liquid interfaces. This is in contrast 
with molecular surfactants, which adsorb to and desorb from the interface 
in a dynamic equilibrium with an excess of surfactant in one of the liquid 
phases. Hence, in the case of emulsification in the presence of colloidal 
particles, coarsening will proceed until the surface area of the droplets can 
only just accommodate all of the particles. The jammed particles then 
halt further phase separation, because the energy barrier preventing their 
expulsion from the liquid-liquid interface is too high. 

Though not extensively, the microstructure of particle-stabilised emul- 
sions has been studied quantitatively, in the case of fabrication by direct 
mixing. Arditty et al. considered dense emulsions of mm-scale droplets, 
formed via limited coalescence and stabilised by clusters of silica nanopar- 
ticles [1| . Note that these samples did not have a surplus of colloids in the 
continuous phase — all particles were interfacial from the start of limited 
coalescence onwards. Using video observations, they found that the aver- 
age droplet diameter rapidly increases initially, after which it saturates at a 
limiting value. The corresponding size distributions are monomodal and sur- 
prisingly narrow, which is attributed to coarsening via coalescence. In short, 
small droplets coalesce fast due to their relatively large surface-to-volume 



ratio (low particle coverage), while large droplets coalesce slowly because 
film thinning, shape relaxation and local rearrangements occur within time 
scales that increase with droplet size. Combined, these effects limit droplet 
polydispersity in their emulsions to 10-20%. 

More recently, Schelero et al. considered direct-mixing emulsions in 
which the micron-sized droplets were stabilised by patches of catanionic 
crystals [12]. Ageing behaviour at room temperature was characterised us- 
ing single-particle light scattering over a period of nine months. Regardless 
of sample composition and age, they found droplet-size distributions with 
three distinctive peaks, which they attribute to 1) pure catanionic crys- 
tals, 2) oil droplets covered with a monolayer of catanionic pairs and 3) oil 
droplets with larger catanionic crystals in addition to the monolayer. Note, 
however, that their emulsion preparation procedure differs significantly from 
that of Arditty et al. [4J. In particular, Schelero et al. added two kinds of 
surfactant that form crystals upon cooling just after mixing all of the in- 
gredients at an elevated temperature. It is not clear to what extent the 
properties of these emulsions are influenced by 1) the presence of surfac- 
tants during emulsification and 2) the lack of solid particles at the start of 
macroscopic phase separation. 

In this paper, we report a quantitative confocal-microscopy study of 
emulsification upon demixing in a binary liquid containing colloidal particles. 
To the best of our knowledge, the microstructure of particle-stabilised emul- 
sions prepared by this route has not yet been characterised quantitatively. 
As emulsification is induced through a controlled temperature quench, the 
process can be imaged in situ and in real time. Obtaining droplet-size dis- 
tributions remains challenging, especially near droplet nucleation, where 
the two demixing liquids are still chemically similar, resulting in an in- 
herently low (fluorescence) imaging contrast. Hence, we have developed an 
image- analysis method that relies on structure factors obtained from discrete 
Fourier transforms of individual frames in confocal image sequences. Radi- 
ally averaging the squared modulus of these Fourier transforms before peak 
fitting allows extraction of dominant length scales over the entire tempera- 
ture range of the quench, even when droplets are not yet discernable in the 
images, i.e. just after nucleation. Thus analysed, our confocal-microscopy 
data implies that our final emulsions have bimodal droplet-size distributions. 

The rest of this paper is organised as follows. Below, in section [2l we 
start with a description of our experimental procedures, notably the image- 
analysis method that we have developed. We continue by presenting our 
main results in section [31 focussing first on the final emulsions and then on 
emulsification. In section [^ we first provide an explanation for the observed 
bimodal droplet-size distributions, after which we discuss emulsion forma- 
tion and stability. Finally, in section [SJ we draw conclusions and suggest 
future work. 



2 Materials and methods 

Sample preparation and characterisation (sections I2.1fl2.3p follow methods 
developed in reference [S]. 

2.1 Sample preparation 

The experiments were performed with silica colloids in the binary mixture 
hexane-methanol [7]. Hydrophobic, fumed silica particles were used as re- 
ceived (Degussa, AEROSIL R812). This aerosil consists of roughly spherical 
primary particles with an average size ~ 7 nm, which are clustered together 
into 0.4 fiui diameter fractal aggregates; these aggregates, as a rule, cannot 
be broken down further |131ll4j. Hexane (Fluka, > 99%), methanol (Fisher 
Scientific, 99.99%) and the fiuorescent dye Nile Red (Sigma) were used as 
received. Nile Red was usually dissolved in the methanol at concentrations 
of 3.7 • 10~^ to 2.9 • 10^^ M. Sample mixtures were typically prepared at a 
hexane/methanol volume ratio of 65/35 (61/39 w/w) and a silica content 
of 2.0 vol% (5.4 wt%). All amounts were determined by weighing. The 
particles were dispersed using an ultrasonic processor at 6 W for 2 minutes 
(Sonics, Vibra-cell), with part of the vial immersed in water, followed by 10 
s of vortex mixing. Evaporation losses during sonication led to a maximum 
deviation of 2.5%-point in the hexane/methanol volume ratio. 

2.2 Sample transfer and cooling 

Sample cells were rectangular and had a 1 mm internal path length (Starna 
Scientific Ltd). Sample mixtures in the single-fiuid phase were transferred 
to these cells in an incubator at approximately 40 °C (Stuart, SI60). They 
were then quickly transferred to a modified hotstage at 40.0 °C (Linkam 
Scientific LTS350). Emulsification was studied while cooling from 40.0 °C 
to 22.0 °C at 5.0 °C • min~^ (figure [D^b)). During an experimental session, 
samples were typically subjected to several heating and cooling cycles. After 
a heating cycle, they were vigorously shaken by hand in the incubator at 
~ 40 °C to ensure proper particle re-dispersion. For samples with less than 
1.0 vol% of silica (figure[2)^c) only), we deviated from the standard procedure 
described above: they were quenched from ~ 40 °C by putting them on a 
metal surface at ~ 22 °C f9]. 

2.3 Sample characterisation by confocal microscopy 

The samples were studied with confocal laser scanning microscopy in re- 
flection, fluorescence and/or transmission. A Nikon ECLIPSE E800/TE300 
upright/inverted microscope was used in conjunction with a BioRad Radi- 
ance 2100 scanning system. The 457 nm line of an Ar-ion laser was used for 
reflection, while a 543 nm HeNe laser was employed to excite the Nile Red. 



Emission filters were used as appropriate. Visual inspection and confocal 
microscopy on samples without particles confirmed that the recorded Nile 
Red fluorescence was mainly coming from the methanol-rich phase. Owing 
to the hotstage, a Nikon Plan Fluor Extra Long Working Distance 20x/0.45 
NA objective with an adjustable correction collar was used. 

2.4 Image analysis 

For quantitative analysis, confocal fluorescence images were enhanced and 
analysed using the IDL software package (RSI v6.3). First, to correct for 
non-uniform brightness, IDL's 'adaptive histogram equalization' function 
was applied to our images. Subsequently, their 2D Fast Fourier Transforms 
were calculated and the modulus of these 2DFFTs was squared. To sim- 
plify the fitting procedure (see below), spatial frequencies close to 0, related 
to the finite size of the images, were removed using a Butterworth high- 
pass filter; within each analysed confocal image sequence, the values of the 
filter parameters were the same for all frames. After that, the |2DFFT| 
arrays were radially averaged and the resulting structure factors S (q) were 
smoothed using a typical box width of 3 to 5 reciprocal pixels. 

Using an iterative Levenberg-Marquardt algorithm, peaks in experimen- 
tal structure factors S (q) were fitted to an 'Extreme function', 

y = y,+^.eH-^- + l), (2) 

where 

z = ^^^ , (3) 

w 

y = S{q), Ho is the offset, A is the peak amplitude, Xc is the peak centre 
and w is its width. Because the images in our confocal time series are 2D 
confocal slices through a 3D emulsion, we define the polydispersity PD2D of 
the population of cross-sectional disks corresponding to a particular peak in 
S{q): 

PD2D = f ^^) • 100% , (4) 

\Xc -Wj 

which roughly corresponds to the relative Half Width at Half Maximum of 
the fitted peaks. The pixel size Aq x Aq of the 2DFFT images was calibrated 
using 

Aq fim'^ = , (5) 

L J iV • Au [jj,m\ 

where Am x Au is the pixel size of the corresponding N x N confocal- 
microscopy image. Finally, we define the dominant length scale in an image 
as 



Peak fitting allowed for batch processing of images in time series of 
demixing by first analysing the last image (A^ — 1), using the resulting pa- 
rameters as an initial parameter estimate for analysing image {N — 2) and 
so on, thus 'backtracking' peaks in S {q) as a function of temperature. While 
this avoids the need to enter trial parameters for each image, it does result 
in spurious fits above the nucleation temperature, i.e. when the sample is in 
the single-fluid phase. 

3 Results 

3.1 Final emulsions 

To facilitate image analysis later on, we first summarise our previous re- 
sults concerning the qualitative structure of hexane-methanol/silica emul- 
sions that had been formed by quenching from the single- fluid phase [S]. 
Figures [2ja) and (b) show the confocal fluorescence and the corresponding 
confocal reflection images of such an emulsion, just after the quench target 
temperature had been reached. The droplets contain the less dense, hexane- 
rich phase [llj . while the continuous phase is rich in methanol. Although 
the silica particles mainly partition into the continuous phase, the confo- 
cal reflection image in figure [2]^c) implies that our emulsions are stabilised 
by interfacially trapped particles, leaving an excess of colloids in the con- 
tinuous phase if the particle volume fraction ^ 0.1%. This implication is 
corroborated by figure [2]^e) , as it seems unlikely that colloids in the contin- 
uous phase alone can stabilise our emulsions for at least 24 minutes at silica 
volume fractions as low as 1.0 vol%. 

To complement these qualitative observations, we calculate structure 
factors S (q) from confocal fluorescence images, thereby revealing that our 
emulsions may well be characterised by bimodal droplet-size distributions. 
Figure EJ^d) shows the structure factor corresponding to figure EJ^a); the 
sample had just reached the quench target temperature. The main peak at 
low q is related to the average distance between large cross-sectional disks 
in the 2D confocal slices of our 3D emulsions; a value of (19.4 it 0.6) /im is 
obtained from the corresponding fit in this case. As the cross-sectional disks 
are (nearly) close-packed, this should roughly correspond to their average 
(apparent) diameter. Note that figure [2]^d) presents a 'worst-case scenario' 
— the width of the main S (q) peak is generally fitted well at half height, but 
the peak asymmetry is usually captured more accurately. Considering the 
relatively high polydispersity PD2D ~ 62% (equation ([3|)), it is unlikely that 
the second peak at higher q stems from the corresponding form factor |15] . 
Instead, we attribute it to the average distance between two small droplets 
or between a small and a large one. In either case, its presence implies a 
bimodal size distribution. The shoulder on the left side of the main peak in 
figure [2{d) is probably an artefact of limited resolution — it may correspond 




0.15. 



0.10- ,jf 



55" 



,A 



0.0 



data 

-fit 



E 



°"7 V ^ 

0.00 -p^ , , 1 I I 1 



32.0. 



31.5- 



31.0- 



O) 30.5 - ^\« 

0) 



30.0. 



•^.•v-'. 



.■^v-v. 



^.^^-A'.^v 



(d) 



0.1 0.2 



0.3 



23.0 23.2 23.4 23.6 23.8 24.0 

(e) sample age / min 



Figure 2: (a,b) Confocal images of a hexane-in-methanol emulsion stabilised 
by 2.0 vol% of silica: (a/b) fluorescence/reflection reveals the methanol-rich 
phase/silica (white). Images recorded at 22.0 °C, following a quench from 
40.0 °C at 5.0 °C-min~^. (c) Confocal reflection image showing interfacially 
trapped particles (0.1 vol% silica), (d) Structure factor S (q) corresponding 
to (a); black arrows point to S (q) peaks, (e) Dominant length scale L 
versus time in an emulsion, containing 1.0 vol% of silica particles, at 22.0 °C 
{AppendixlQj . Images were recorded at 14 fiui from the top or (c) 0.79 mm 
from the bottom of the sample. Scale bars: 100 /um. Panel (c) reproduced 
by permission of The Royal Society of Chemistry [S] . 



to pairs of neighbouring droplets that appear as dumbbells to the software. 
Having considered the structure of the final emulsions, we will next use this 
structure-factor analysis to explore the emulsification process. 

3.2 Emulsification 

Confocal fluorescence image sequences of the emulsification process provide 
a clue as to the origin of the observed bimodal droplet-size distributions. 
Figure [3] shows a snapshot montage extracted from such a sequence {Ap- 
pendix\Q. At the start of the quench (figure EH^a)), the fluorescence signal is 
flat across the frame, which is consistent with a binary-liquid/particle sys- 
tem in its single-fluid phase. At the nucleation temperature (figure O^b)), 
as determined from an intensity drop in the transmission channel, the flu- 
orescence image is still flat, probably because 1) our confocal resolution is 
not sufficient to resolve the nucleated droplets [16j and 2) the two demixing 
liquids still have similar chemical compositions, limiting the (fluorescence) 
contrast. Once droplets have become discernable (flgures Ws) and (d)), 
phase separation appears to have yielded a monomodal size distribution. 
Upon further cooling, however, larger droplets start popping up in between 
consecutive confocal frames (figures [3)^f)-(i)). We attribute this to coarsen- 
ing, at this stage, being dominated by coalescence together with creaming, 
clearly resulting in a multimodal size distribution. 

To quantify this coarsening behaviour, experimental structure factors 
were fitted to a functional form (equations ([2]) and ([3])), thereby allowing 
batch processing of confocal images. Figure U)^ a) shows a graph of the domi- 
nant length scale L as a function of temperature (equation dG])), the features 
of which can be explained by comparing it to the corresponding image se- 
quence in figure [3l Before droplet nucleation ( A-B) , the confocal images 
contain no dominant length scale, apart from the image size and the pixel 
size. At the nucleation temperature (B), determined independently from a 
fall in transmission intensity, our software picks up a reasonable length scale, 
even though no droplets are discernable in the corresponding confocal image 
(figure [3)^b))! Just after droplet nucleation (B-D), L increases quickly as the 
sample is quenched deeper into the two-phase part of the phase diagram. 
From (D-E), the growth rate falls, but it picks up again at temperatures be- 
low (E)jj which we ascribe to coalescence taking over as the main coarsening 
mechanism (figure [3|^e-i) ) . 

As indicated by the vertical line at 22.0 °C in the graph in figureSJa), our 
emulsions do not necessarily stop coarsening immediately after the quench. 
This is revealed more clearly in figure Hl^b), in which the same data has 
been re-plotted as a function of time. For this sample, just after the target 



^In similar graphs for emulsions containing 2.5 and 3.5 vol% of silica (not shown here), 
the increase in the rate of droplet growth at temperatures below (E) is less pronounced 
(or absent). 
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Figure 3: Selected xy-frames, analysed in figure[3]and centrally cropped here 
(512 X 512 — )• 256 x 256 pixels), from a confocal fluorescence image sequence 
of emulsification in a hexane-methanol mixture containing 2.0 vol% of silica 
particles (see Appendix fCl for movie and Appendix W\ for structure factors); 
panel (h) is a digital zoom of figure [2||a) . Starting at t = s, the sample 
was quenched from 40.0 °C to 22.0 °C at 5.0 °C • min~ , resulting in the 
nucleation and growth of hexane-rich droplets in a methanol-rich continuous 
phase. Images were recorded at 14 /im from the top of the sample. Scale 
bars: 100 /im. 
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Figure 4: Probing the dominant length scale in the image sequence shown 
in figure [3l the uppercase letters 'A' to T correspond to the lowercase panel 
labels in figure [3l Note that the images in figure [3] have been centrally 
cropped for clarity (512 x 512 — )• 256 x 256 pixels), but the original images 
have been analysed for the graphs above. (a,b) Dominant length scale and 
(c) polydispersity PD2D as a function of (a,c) temperature and (b) time. 
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temperature was reached, the coarsening rate started to decrease. However, 
it took about a minute before the dominant length scale reached a plateau 
value. Judging from experiments with a thermocouple embedded in a sam- 
ple without colloids, the time lag between sample and hotstage temperatures 
during a 5.0 °C • min~ quench may be as large as 50 s in our set-up, which 
could explain the residual coarsening. However, in similar graphs for emul- 
sions containing 2.5 and 3.5 vol% of silica (not shown here), the dominant 
length scale does not reach a plateau value within the same time interval, so 
it is unclear whether our emulsions are stable immediately after the quench. 
To probe the stability of our emulsions on a longer time scale, the dominant 
length scale L in a 23-minute old emulsion was extracted from a one-minute 
time series at constant temperature (figure [2]||e)). This graph demonstrates 
that our emulsions are stable on experimental time scales, for L is virtu- 
ally constant; the two jumps near t = 23.2 and 23.8 min correspond to 
re-arrangements coinciding with coalescence events {Appendices\M and K^ . 

4 Discussion 

4.1 Emulsion structure 

The two peaks in the structure factor in figure [2](d) imply that our emulsi- 
fication route leads to bimodal particle-stabilised emulsions, which we have 
attributed to coalescence together with creaming (figure [3]) . The underlying 
mechanism we propose is schematically depicted in figure Early in the 
emulsification process, the liquid-liquid interface is not fully coated with par- 
ticles (figure [SK a)). Though there is a surplus of colloids in the continuous 
phase (figure [2](b)), the droplets are mainly stabilised by interfacial parti- 
cles (figure Wis) ) ; so a partial coating is not sufficient to prevent them from 
merging. As volume is conserved upon coalescence, the surface-to-volume 
ratio V decreases when two spherical droplets of radii Ri and R2 merge to 
form a spherical droplet of radius i?3 (figure [5l|a) — T-lSKb)): 



T^l 



1^3 



47rfif _ _3_ 

(4/3)7rRf Ri 



(A) 



1+ ^ 



vv^^v^y J 



(7) 



Since the number of interfacial particles is (approximately) conserved during 
coalescence, their area fraction increases, which yields droplets with more 
complete particle coatings, i.e. a more stable emulsion (figure[5^b)). As the 
quench progresses, unmerged droplets can sweep up more particles from the 
surplus in the continuous phase, thereby completing their coatings (figure 
[5]^c)). Note that this last step does not occur during limited coalescence, 
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Figure 5: 2D schematic explaining biniodal droplet-size distributions in 
emulsions formed upon quenching a binary-liquid/particle mixture from its 
single-fluid phase, (a) Just after droplet nucleation, the droplets (black) are 
only partly covered with particles (grey), (b) When two droplets merge, the 
resulting droplet has a smaller surface-to- volume ratio (equation dZ])), i.e. it 
has a denser particle coating and is thus more stable, (c) Droplets that have 
not merged can sweep up additional particles from the continuous phase 
(white), thus stabilising themselves vs coalescence. This emulsification pro- 
cess results in two 'generations' of stable droplets with different average 



sizes. 



for all particles are interfacial in that process. This may very well be why 
Arditty et al. reported monodisperse emulsions [Ij, whereas we observe two 
'generations' of droplets with different average diameters. 

At this point, it is worthwhile recalling that our confocal images are 
2D slices at a distance d = 14 //m from the top of the sample and, hence, 
they do not necessarily represent the whole 3D emulsion. However, we will 
argue below that, in our case, a bimodal size distribution for the 2D cross- 
sectional disks implies a bimodal size distribution for the corresponding 3D 
droplets, even though the exact shapes of their graphs may differ. First of 
all, our emulsions appear to be random packings of polydisperse spherical 
droplets in 3D [9]. It is unlikely that a 2D slice through such a 3D emulsion 
would yield a bimodal disk-size distribution from a monomodal droplet-size 
distribution. Secondly, assuming that all droplets in the imaging plane have 
a 3D radius R > d/2 and are touching the sample-cell wall, which is not 
unreasonable during the late stages of emulsification in which a compact 
cream has formed at the top of the sample, one can write down the relation 
between R and the 2D disk radius r: 



R 



2d 



(8) 



As r > 0, equation ^ expresses a one-to-one correspondence between r 
and R, i.e. if the 2D disks are characterised by a bimodal size distribution, 
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then so are the 3D droplets. Finally, 3D simulations on nucleation and 
growth of droplets in similar systems have also found bimodal droplet-size 
distributions which can be ascribed to coalescence [T7|; our experimental 
results seem to confirm these simulations. 

4.2 Emulsification 

Our emulsification process begins with the nucleation of hexane-rich droplets 
upon cooling the sample from the single- fluid phase into the two-phase region 
of the phase diagram at an off-critical composition (figures [3]^b/c) and[Hb)). 
Intriguingly, nucleation may happen at temperatures above the binary-liquid 
binodal, at 36.5 °C instead of 33.2 °C in this case (figures[H^a) and[Hb)) [7J, 
though the actual value of the nucleation temperature varies from sample 
to sample. This shift is partly due to a temperature lag between sample 
and hotstage, which can be as large as 0.5 °C in our set-up. In addition, 
contamination with only a few vol% of water, which is easily picked up from 
the air, increases the consolute temperature of the binary mixture by several 
°C [18] o Finally, colloidal particles may raise the binodal of a binary liquid, 
by as much as 1.0 °C/vol% in this system [9j. 

An elevated nucleation temperature that varies with particle volume frac- 
tion is likely to imply heterogeneous nucleation [19] . Indeed, it is essentially 
impossible to prevent heterogenous nucleation when quenching off-critical 
mixtures [16j. However, given the strong preference of our silica particles 
for the continuous methanol-rich phase, it seems unlikely that the hexane- 
rich phase would nucleate onto these colloids to form droplets. Instead, 
we propose that the solid particles acquire a methanol-rich layer due to 
adsorption/wetting, thereby partially depleting the volume in between the 
particles of methanol, i.e. locally increasing hexane supersaturation. This 
would indeed result in particle-promoted nucleation of hexane-rich droplets 
at temperatures above the binodal of the binary liquid itself. 

During the initial stage of emulsification (B-D), liquid-liquid demixing 
appears to yield a monomodal droplet-size distribution. Though merely 
speculation, this could be related to the inherently narrow size distribution 
resulting from free growth of droplets in a dilute emulsion with a continu- 
ous phase at quasi-constant supersaturation ( Appendix W\i. Indeed, figure 
[3{c/d) shows that our emulsions are initially quite dilute. The correspond- 
ing volume fraction of droplets (</>dropict) can be estimated using the lever 
rule, 



ooSF _ MR 
^M ^M 



^T - ^T 



(9) 



'/'droplet 

where (/?|j is the volume fraction of methanol in the single-fluid phase and 



^Contamination with water can be minimised by storing the hexane over molecular 
sieves and the methanol under dry nitrogen gas until the day of the experiment. 
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ify^ is the volume fraction of methanol in the hexane(methanol)-rich 

phase [2U]. At the same temperature difference from the binodal as in figure 
[3^d), equation ([9]) predicts i;^dropiet ~ 0.32, which is well below the onset of 
crystallisation for monodisperse hard spheres |21j . In addition, we expect 
the methanol-rich continuous phase to be constantly supersaturated with 
hexane, as our quench is continuous. However, given a 1.20 //m x 1.20 yum 
pixel size, a critical nucleation radius ~ 0.03 ^m [161 ES] and a particle- 
cluster size ~ 0.4 fiui (section 12. ip . it is unlikely that we have probed the 
free-growth regime itself — we may have observed a remnant thereof. 

During the intermediate stage of emulsification (D-E) , droplet growth in 
the imaging plane, which is near the top of the sample, slows down (figure 
[Ula)). This may be explained by the disappearance of (three) buoyancy- 
related contributions to droplet growth as the cream develops. Firstly, part 
of the growth in the initial stage may have been caused by large droplets 
rising into the imaging plane by steady creaming [23j, an effect that will 
contribute less and less as the observed droplet layer densities. Secondly, 
as the cream becomes more compact, the local volume fraction of methanol 
decreases, i.e. the system moves towards the hexane-rich branch of the phase 
diagram (figure mb)) [9]. This means that, in/near the imaging plane, the 
hexane supersaturation of the continuous phase decreases, which slows down 
droplet growth. Finally, creaming will induce convective flows that may 
contribute to droplet growth at least as much as diffusion does [111 [T6] . 

During the final stages of emulsification, at temperatures below (E), 
coalescence takes over as the main coarsening mechanism (figure [3] and Ap- 
pendices\^ and\Q) . The reasons for this are twofold. First of all, buoyancy 
has resulted in a compact cream at the top of the sample, in which many 
droplets are now close enough to merge. Secondly, as the temperature has 
now significantly fallen below the binodal, the interfacial tension between 
the two liquid phases has increased substantially [H]. This drives emulsion 
coarsening through coalescence as the system attempts to reduce the total 
free-energy cost associated with liquid-liquid contact area. In addition, in- 
creasing interfacial tension promotes particle attachment (equation ([I|)). As 
explained in section 14. H the presence of colloidal particles on the surface 
of the droplets eventually stabilises them, resulting in bimodal droplet-size 
distributions (figure [5]). 

4.3 Emulsion stability 

Particle-stabilised emulsions can be long lived [T]; so why do our emulsions 
slowly coarsen after the quench has reached its target temperature (figures 
Hl^b) and[2]^e))? Our system is comprised of droplets stabilised by interfacial 
particles (figure [2]^c)). These fumed silica colloids are fractal clusters of 
~ 7 nm primary particles and have a total diameter ~ 0.4 //m |13[ I14j . 
As equation ([1]) predicts a maximum free energy of detachment AG^ ~ 
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10"* ^B^l-oom for 0.4 fiui particles [UITT], it seems unlikely that they could be 
expelled from the liquid-liquid interface. However, the surplus of particles in 
the continuous phase suggests that the particle-liquid-liquid wetting angle 
9 is some way away from 90° (figure EJb)) [9]. Moreover, as the fumed 
silica particles are fractal clusters [13^ I14j . they cover less interfacial area 
than equally sized spheres do and hence experience a shallower energy well 
than that described by equation ([T]). Finally, for partially miscible liquids, 
Ostwald ripening may also be quite severe [21]. Taken together, these effects 
may explain the observed slow coarsening. 

5 Conclusions 

In this paper, we have presented a quantitative confocal-microscopy study of 
emulsification upon demixing in a binary liquid containing colloidal particles. 
Our analysis method relies on structure factors S {q) calculated from individ- 
ual images in confocal fluorescence time series via discrete 2D Fast Fourier 
Transforms. Radially averaging the squared modulus of these 2DFFTs be- 
fore peak fitting allows extraction of droplet size and variance over the entire 
temperature range of the emulsifying quench. 

First of all, our investigation implies that our particle-stabilised emul- 
sions are characterised by bimodal droplet-size distributions |j We attribute 
this to coalescence together with creaming being the main coarsening mech- 
anism during the late stages of emulsification and we support this claim with 
(direct) confocal-microscopy observations. Secondly, we have quantitatively 
characterised the emulsification process itself. Our results imply that it 
starts with particle-promoted nucleation, possibly followed by free growth 
of droplets. During the intermediate stage, in which we observe the forma- 
tion of a compact cream near the top of the sample, coarsening slows down. 
This may be explained by the disappearance of (three) buoyancy-related 
contributions to droplet growth: 

1. large droplets no longer move into the imaging plane from below; 

2. the increasing packing fraction of hexane-rich droplets in the creamed 
emulsion reduces hexane supersaturation of the methanol-rich contin- 
uous phase; 

3. fading of convective flows, which may contribute to droplet growth at 
least as much as diffusion does. 

During the late stages of emulsification, at temperatures significantly below 
the binodal, the close packing and the increased interfacial tension result in 
emulsion coarsening through coalescence. We have shown that the presence 



•^Preliminary experiments suggest that slower quenches at 0.1 °C ■ min ^ result in 
multimodal rather than bimodal particle-stabilised emulsions. 
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of interfacial particles eventually stabilises the emulsions on experimental 
time scales (~ minutes). However, fumed silica tends to be organised in 
fractal clusters, which are not efficient in covering liquid-liquid interface. 
This may explain the slow coarsening of our emulsions at constant temper- 
ature. 

To conclude, phase separation is a promising route towards the controlled 
formation of particle-stabilised emulsions. The experimental conditions dur- 
ing emulsification can be controllably varied in-situ and in real time. A more 
detailed study, involving the tracking of individual droplets in microscopy 
images, should lead to a better understanding of the relation between quench 
conditions and emulsion structure. As such, binary liquids containing col- 
loidal particles are an ideal model system for studying emulsification and 
could potentially allow for tailor-made particle-stabilised emulsions. 
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A Imaging coalescence 
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Figure 6: Rows: confocal fluorescence image sequences of coalescence events 
in hexane- methanol emulsions, (a-f) Starting at t = s, the sample was 
quenched from 40.0 °C to 22.0 °C at 5.0 °C • min^^. These images were 
taken from the same series as those in figure 3 of the main text, (g-i) The 
quench target temperature was reached at t = s, after which the sample 
was held at 22.0 °C. These images were taken from the series corresponding 
to figure 2(e) of the main text (Appendix l^ . Images were recorded at 14 
/im from the top of the sample. Scale bars: 100 /xm. 
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B Monodisperse emulsions from free growth of 
droplets 

To understand why monodisperse, and hence monomodal, emulsions are 
an inherent feature of droplet nucleation followed by a free-growth regime 
[151 [IS]j we turn to the growth equation for an isolated droplet in a super- 
saturated continuous phase 



dp 

dr 


= P-' (l - P-' 




P =W 

*diff 



(10) 

in which 

„ — JL 

(11) 

Here, R is the droplet radius, t the time after nucleation and tdiflf a diffusion 
time scale. -Re is the critical radius: nucleated droplets with R < R^ will 
shrink and disappear, while those with R > R^ will grow. For R ^ Re, 
equation PU]) reduces to 

P=V2^, (12) 

whatever the initial distribution of droplets [13 dH]- In words: droplets 
smaller than Re disappear in order to reduce the free energy of the system, 
while droplets much larger than Re have locally depleted the continuous 
phase of the nucleating phase, thereby reducing supersaturation and slowing 
their growth. 

Note that it is not unreasonable to assume that equations (fTO]l and (fTTj) 
apply in our case as 1) our emulsions are initially quite dilute, 2) the con- 
tinuous phase is supersaturated (before creaming) and 3) the droplets may 
very well nucleate in between rather than on the colloidal particles. Even in 
the case of heterogeneous nucleation, the above argument may hold |15[I16). 
as the critical radius Re is still a well-defined concept. If one assumes that 
the embryo nucleating onto a spherical seed particle is itself a portion of a 
different sphere of radius R, then Re is the critical radius of this embryonic 
sphere. Notably, the value of Re is the same as in the case of homogeneous 
nucleation, even though the volume of the embryo will differ from that of a 
homogenously nucleated droplet, as all parts of the embryo surface must be 
in equilibrium with the metastable continuous phase 
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C Movie guide 

See http://tinyurl.com/382tmz2 

• Thijssen_2010-EQA_SI-Movie2e-2xrt.avi: corresponds to the graph in 
figure 2(e) of the main text. Hexane-in- methanol emulsion stabilised 
by 1.0 vol% of silica particles. The fluorescence reveals the methanol- 
rich phase (white). At the start of this movie, the sample had already 
been held at a constant temperature of 22.0 °C for 23 minutes. Images 
were recorded at a depth of 14 //m from the top of the sample. Image 
dimensions: 307.97 ^m x 307.97 ^m. The movie is at 2 x real time; it 
has been resized from 512 x 512 to 384 x 384 pixels and subsequently 
JPEG compressed to reduce file size. 

• Thijssen_2010-EQA_SI-Movie3-6xrt.avi: corresponds to figure 3 of the 
main text. Confocal fluorescence image sequence of emulsification in a 
hexane-methanol mixture containing 2.0 vol% of silica particles. Start- 
ing at i = s, the sample was quenched from 40.0 °C to 22.0 °C at 
5.0 °C • min~^, resulting in the nucleation and growth of hexane-rich 
droplets in a methanol-rich continuous phase. Images were recorded 
at a depth of 14 /im from the top of the sample and at intervals of 3 s. 
Image dimensions: 307.97 //m x 307.97 /im. The movie is at 6x real 
time and has been JPEG compressed to reduce file size; reproduced 
by permission of The Royal Society of Chemistry [9]. 
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D Structure factors for image sequence in figure 
3 (main text) 
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Figure 7: Structure factors S (q) corresponding to figure 3(a-i) of the main 
text, which shows a confocal fluorescence image sequence of eniulsification in 
a hexane-methanol mixture containing 2.0 vol% of silica particles {Appendix 
\Q). The emergence and growth of the main(secondary) peak reflects the 
(subsequent) emergence of a dominant (secondary) length scale. Starting at 
t = s, the sample was quenched from 40.0 °C to 22.0 °C at 5.0 °C • min^^, 
resulting in the nucleation and growth of hexane-rich droplets in a methanol- 
rich continuous phase. The analysed images were recorded at 14 /xm from 
the top of the sample. 
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